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Abstract 

In this paper we propose the time-dependent Hamiltonian form of human biomechanics, as 
a sequel to our previous work in time-dependent Lagrangian biomechanics [T]. Starting with 
the Covariant Force Law, we first develop autonomous Hamiltonian biomechanics. Then we 
extend it using a powerful geometrical machinery consisting of fibre bundles and jet manifolds 
associated to the biomechanical configuration manifold. We derive time-dependent, dissipa- 
tive, Hamiltonian equations and the fitness evolution equation for the general time-dependent 
human biomechanical system. 

Keywords: Human biomechanics, covariant force law, configuration manifold, jet manifolds, 
time-dependent Hamiltonian dynamics 

1 Introduction 

Most of dynamics in contemporary human biomechanics is autonomous (see [2j [3J SI [5j El [3 [HI [Si 
1101 111] ). This approach works fine for most individual movement simulations and predictions, in 
which the total human energy dissipations are insignificant. However, if we analyze a lOOm-dash 
sprinting motion, which is in case of top athletes finished under 10 s, we can recognize a significant 
slow-down after about 70 m in all athletes - despite of their strong intention to finish and win the 
race, which is an obvious sign of the total energy dissipation. This can be seen, for example, in 
a current record-braking speed-distance curve of Usain Bolt, the world-record holder with 9.69 s, 
or in a former record-braking speed-distance curve of Carl Lewis, the former world-record holder 
(and 9 time Olympic gold medalist) with 9.86 s (see Figure 3.7 in [H]). In other words, the total 
mechanical energy of a sprinter cannot be conserved even for 10 s. So, if we want to develop a 
realistic model of intensive human motion that is longer than 7-8 s, we necessarily need to use the 
more advanced formalism of time-dependent mechanics. 

In this paper, as a sequel to our previous work in time-dependent Lagrangian biomechanics [T], 
we use the covariant force law in conjunction with the modern geometric formalism of jet manifolds 
to develop general Hamiltonian approach to time-dependent human biomechanics in which total 
mechanical energy is not conserved. 
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2 The Covariant Force Law 



Autonomous Hamiltonian biomechanics (as well as autonomous Lagrangian biomechanics) , based 
on the postulate of conservation of the total mechanical energy, can be derived from the covariant 
force law [3J [3J [H [5] , which in 'plain English' states: 

Force 1-form = Mass distribution x Acceleration vector-field, 

and formally reads (using Einstein's summation convention over repeated indices): 

F l = m ij a j . (1) 

Here, the force 1-form F{ — Fj(t, q,p) = F((t, q, q), (i — 1, n) denotes any type of torques and 
forces acting on a human skeleton, including excitation and contraction dynamics of muscular- 
actuators [TH [13j [12] and rotational dynamics of hybrid robot actuators, as well as (nonlinear) 
dissipative joint torques and forces and external stochastic perturbation torques and forces [BJ. 
rriij is the material (mass-inertia) metric tensor, which gives the total mass distribution of the 
human body, by including all segmental masses and their individual inertia tensors. a J is the total 
acceleration vector-field, including all segmental vector-fields, defined as the absolute (Bianchi) 
derivative v % of all the segmental angular and linear velocities v % = x 1 , (i = 1, n), where n is 
the total number of active degrees of freedom (DOF) with local coordinates (x l ). 




Figure 1: Biomechanical configuration manifold defined as an anthropomorphic product tree con- 
sisting of Euclidean motion groups SE(3) (described in the text below). 

More formally, this central Law of biomechanics represents the covariant force functor 
defined by the commutative diagram: 
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TT*M 
Fi=p. 



T 



■TTM 



a = v 

T*M={x i ,p i } TM = {x\v i } 
M = {x { } 



(2) 



Here, M = M n — {x 1 . (i = 1, ...,n)} is the biomechanical configuration n— manifold (Figure 
[T]), that is the set of all active DOF of the biomechanical system under consideration (in general, 
human skeleton), with local coordinates (x l ). 

The right-hand branch of the fundamental covariant force functor T„ : TT*M — ► TTM de- 
picted in @ is Lagrangian dynamics with its Riemannian geometry. To each n— dimensional (nD) 
smooth manifold M there is associated its 2nD velocity phase-space manifold, denoted by TM 
and called the tangent bundle of M. The original configuration manifold M is called the base of 
TM. There is an onto map tt : TM — ► M, called the projection. Above each point x G M there 
is a tangent space T X M = tt (x) to M at x, which is called a fibre. The fibre T X M C TM is 
the subset of TM, such that the total tangent bundle, TM = j | T X M, is a disjoint union of 

m£M 

tangent spaces T X M to M for all points x G M. From dynamical perspective, the most important 
quantity in the tangent bundle concept is the smooth map v : M — > TM, which is an inverse to 
the projection tt, i.e, tt o v = Wm, tt(v(x)) — x. It is called the velocity vector-field v l — x 1 Its 
graph (x,v(x)) represents the cross-section of the tangent bundle TM. Velocity vector-fields are 
cross-sections of the tangent bundle. Biomechanical Lagrangian (that is, kinetic minus potential 
energy) is a natural energy function on the tangent bundle TM. The tangent bundle is itself a 
smooth manifold. It has its own tangent bundle, TTM. Cross-sections of the second tangent 
bundle TTM are the acceleration vector-fields. 

The left-hand branch of the fundamental covariant force functor T„ : TT*M — > TTM depicted 
in ^ is Hamiltonian dynamics with its symplectic geometry. It takes place in the cotangent bundle 
T*M ro b, defined as follows. A dual notion to the tangent space T X M to a smooth manifold M at a 
point x = (x l ) with local is its cotangent space T*M at the same point x. Similarly to the tangent 
bundle TM, for any smooth nD manifold M, there is associated its 2nD momentum phase-space 
manifold, denoted by T* M and called the cotangent bundle. T* M is the disjoint union of all its 
cotangent spaces T*M at all points x G M, i.e., T*M = T*M. Therefore, the cotangent 

xeM 

bundle of an n— manifold M is the vector bundle T*M = (TM)*, the (real) dual of the tangent 
bundle TM. Momentum 1 -forms (or, covector-fields) pi are cross-sections of the cotangent bundle. 
Biomechanical Hamiltonian (that is, kinetic plus potential energy) is a natural energy function on 
the cotangent bundle. The cotangent bundle T*M is itself a smooth manifold. It has its own 
tangent bundle, TT*M . Cross-sections of the mixed-second bundle TT*M are the force 1-forms 
Fi =pi. 

lr This explains the dynamical term velocity phase-space, given to the tangent bundle TM of the manifold M. 
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There is a unique smooth map from the right-hand branch to the left-hand branch of the 
diagram |2]): 

TM 3 {x\v l ) i — ► (x^p 1 ) S T*M. 

It is called the Legendre transformation, or /ifeer derivative (for details see, e.g. [H[S]). 

The fundamental covariant force functor : TT*M — > TTM states that the force 1-form 
Fi = pi, defined on the mixed tangent-cotangent bundle TT*M, causes the acceleration vector- 
field a 1 — v l , defined on the second tangent bundle TTM of the configuration manifold M . The 
corresponding contravariant acceleration functor is defined as its inverse map, T* : TTM — > 
TT*M. 

Representation of human motion is rigorously defined in terms of Euclidean 5 , £ , (3)-groupfH of 
full rigid-body motion in all main human joints 7}. The configuration manifold M for human 
musculo-skeletal dynamics is defined as a Cartesian product of all included constrained SE(3) 
groups, M — Y\j SE(3y where j labels the active joints. The configuration manifold M is coordi- 
nated by local joint coordinates x l (t), i — 1, n = total number of active DOF. The corresponding 
joint velocities x l {t) live in the velocity phase space TM, which is the tangent bundle of the con- 
figuration manifold M. 

The velocity phase-space TM has the Riemannian geometry with the local metric form: 

(g) = ds 2 = gijdx 1 dx\ (3) 

where gij = g t j (m, x) is the material metric tensor defined by the biomechanical system's mass- 
inertia matrix and dx l are differentials of the local joint coordinates x l on M. Besides giving the 
local distances between the points on the manifold M, the Riemannian metric form (g) defines the 
system's kinetic energy: 

T = —gtjX x? , 

giving the Lagrangian equations of the conservative skeleton motion with kinetic-minus-potential 
energy Lagrangian L = T — V , with the corresponding geodesic form 

j t L ± ,-L x , = Q, or ■ Y\,yy I). (4) 

where subscripts denote partial derivatives, while T^ k are the Christoffel symbols of the affine 
Levi-Civita connection of the biomechanical manifold M. 

2 Briefly, the Euclidean SE(3)-group is denned as a semidirect (noncommutative) product (denoted by >) of 3D 
rotations and 3D translations: SE(3) := SO(3) > R 3 . Its most important subgroups are the following: 



Subgroup 




Definition 




50(3), group of rotations 
in 3D (a spherical joint) 


Set of all proper orthogonal 
3 X 3 — rotational matrices 


SE(2), special Euclidean group 
in 2D (all planar motions) 


St 


it of all 3 X 3 — matrici 
cos 8 sin 8 r x 
— sin 8 cos 8 r y 
1 


:s: 


SO(2), group of rotations in 2D 
subgroup of S l £(2)-group 
(a revolute joint) 


Set of all proper orthogonal 
2 X 2 — rotational matrices 
included in SE(2) — group 


R a , group of translations in 3D 
(all spatial displacements) 


Euclidean 3D vector space 
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The corresponding momentum phase-space P = T*M provides a natural symplectic struc- 
ture that can be defined as follows. As the biomechanical configuration space M is a smooth 
n— manifold, we can pick local coordinates {dx 1 , dx 11 } 6 M. Then {dx 1 , dx n } defines a basis 
of the cotangent space T*M, and by writing 9 £ T*M as 9 = Pidx 1 , we get local coordinates 
{x 1 , x n ,pi, on T*M. We can now define the canonical symplectic form wonP = T*M 

as: 

u) = dpi A dx 1 , 

where 'A' denotes the wedge or exterior product of exterior differential formsH This 2— form u> is 
obviously independent of the choice of coordinates {x 1 , ...,x n } and independent of the base point 
{x 1 , ...,x n ,pi, —,p n } S T*M. Therefore, it is locally constant, and so dtu = 00 

If (P, ijj) is a 2nD symplectic manifold then about each point ieP there are local coordinates 
{x 1 , ...,x n ,pi, ...,pn} such that u> = dpiAdx 1 . These coordinates are called canonical or symplectic. 
By the Darboux theorem, u> is constant in this local chart, i.e., du> = 0. 

3 Autonomous Hamiltonian Biomechanics 

We develop autonomous Hamiltonian biomechanics on the configuration biomechanical manifold 
M in three steps, following the standard symplectic geometry prescription (see [H [H O [8] ) : 

Step A Find a symplectic momentum phase-space (P,uj). 

Recall that a symplectic structure on a smooth manifold M is a nondegenerate closecH 2— form 

3 Recall that an exterior differential form a of order p (or, a p—form a) on a base manifold X is a section of the 
p 

exterior product bundle AT*X — > X. It has the following expression in local coordinates on X 

a = ax 1 ...\ p dx Xl A ■ ■ ■ A dx x ? (such that \a\=p), 

where summation is performed over all ordered collections (Ai, X p ). ii p (X) is the vector space of p— forms on a 
biomechanical manifold X. In particular, the 1— forms are called the Pfaffian forms. 

4 The canonical 1— form 6 on T* M is the unique 1— form with the property that, for any 1— form /3 which is a 
section of T* M we have f3*8 = 8. 

Let / : M — » M be a diffeomorphism. Then T* f preserves the canonical 1— form 9 on T*M, i.e., (T* f)*0 = 0. 
Thus T* f is symplectic diffeomorphism. 

5 A p—form (3 on a smooth manifold M is called closed if its exterior derivative d = didx 1 is equal to zero, 

d/3 = 0. 

From this condition one can see that the closed form (the kernel of the exterior derivative operator d) is con- 
served quantity. Therefore, closed p— forms possess certain invariant properties, physically corresponding to the 
conservation laws. 

Also, a p—form j3 that is an exterior derivative of some (p — 1)— form a, 

(3 = da, 

is called exact (the image of the exterior derivative operator d). By Poincare lemma, exact forms prove to be closed 
automatically, 

d/3 = d(da) = 0. 

Since d? = 0, every exact form is closed. The converse is only partially true, by Poincare lemma: every closed 
form is locally exact. 

Technically, this means that given a closed p—form a £ Q P (U), defined on an open set U of a smooth manifold 
M any point m £ U has a neighborhood on which there exists a (p — 1)— form j3 S Cl p ~ 1 (U) such that dj3 = a\u ■ 
In particular, there is a Poincare lemma for contractible manifolds: Any closed form on a smoothly contractible 
manifold is exact. 
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lo on M, i.e., for each x £ M, is nondegenerate, and = 0. 

Let T*M be a cotangent space to M at m. The cotangent bundle T*M represents a union 
U me MT*M, together with the standard topology on T*M and a natural smooth manifold structure, 
the dimension of which is twice the dimension of M. A 1— form 8 on M represents a section 
9 : M -> T*M of the cotangent bundle T*M. 

P = T*M is our momentum phase-space. On P there is a nondegenerate symplectic 2— form 
u> is defined in local joint coordinates x l ,pi S P, U open in P, as co = dx % A dpi. In that case the 
coordinates x l ,pi € U are called canonical. In a usual procedure the canonical 1— form 9 is first 
defined as 9 — p{dx % , and then the canonical 2-form lo is defined as lo = —d9. 

A symplectic phase-space manifold is a pair (P, lo). 

Step B Find a Hamiltonian vector-field Xh on (P,lo). 

Let (P, be a symplectic manifold. A vector-field X : P — > TP is called Hamiltonian if 
there is a smooth function F : P — > R such that = <iP denotes the interior product or 

contraction of the vector-field X and the 2-form lo). X is locally Hamiltonian if ixw is closed. 

Let the smooth real-valued Hamiltonian function H : P — > R, representing the total biome- 
chanical energy H(x,p) — T(p) + V(x) (T and V denote kinetic and potential energy of the 
system, respectively), be given in local canonical coordinates x z ,pi £ U, U open in P. The Hamil- 
tonian vector-field Xjj, condition by ixn 1 - = dH, is actually defined via symplectic matrix J, in 
a local chart U, as 

X H = JWH = (d Pi H,-d x iH), J =(°i J)' ( 5 ) 
where / denotes the n x n identity matrix and V is the gradient operator. 

Step C Find a Hamiltonian phase-flow <j) t of Xh- 

Let (P, ui) be a symplectic phase-space manifold and Xh — JVP a Hamiltonian vector-field 
corresponding to a smooth real- valued Hamiltonian function H : P — > R, on it. If a unique one- 
parameter group of diffeomorphisms 4> t : P — > P exists so that 4f\t=o4> t x — JVH(x), it is called 
the Hamiltonian phase-flow. 

A smooth curve t i— > Pi(t)) on (P, w) represents an integral curve of the Hamiltonian 

vector-field = JVPT, if in the local canonical coordinates 6 P, U open in P, Hamiltonian 
canonical equations hold: 

q l = d Pt H, Pi = -d q iH. (6) 

An integral curve is said to be maximal if it is not a restriction of an integral curve defined on 
a larger interval of R. It follows from the standard theorem on the existence and uniqueness of the 
solution of a system of ODEs with smooth r.h.s, that if the manifold (P, to) is Hausdorff, then for 
any point x — (x l ,pi) G U, U open in P, there exists a maximal integral curve of Xh = JVP, 
passing for t = 0, through point x. In case is complete, i.e., Xh is C p and (P,co) is compact, 
the maximal integral curve of Xh is the Hamiltonian phase-flow <f> t : U — ► P. 

The phase-flow t is symplectic if a; is constant along </> t , i.e., 0J"cj = 
denotes the pu/P6ac4fl of to by </> t ), 

6 Given a map / : X — > X' between the two manifolds, the pullback on X of a form a on X' by / is denoted by 
fa. The pullback satisfies the relations 

f(aAfl = faA/*ft df*a = f*(da), 



iff £ Xh u = 

(£,Xh^> denotes the Lie derivativ^\ of u> upon Xh)- 

Symplectic phase-flow <fi t consists of canonical transformations on (P, u>), i.e., diffeomorphisms 
in canonical coordinates x l ,pi £U, U open on all (P,lj) which leave w invariant. In this case the 
Liouville theorem is valid: <j> t preserves the phase volume on (P, u>). Also, the system's total energy 
H is conserved along <j> t , i.e., H o <f) t = </> t . 

Recall that the Riemannian metrics g — <, > on the configuration manifold M is a positive- 
definite quadratic form g : TM — ► M, in local coordinates a; 1 G U, U open in M, given by © above. 
Given the metrics gij, the system's Hamiltonian function represents a momentum p-dependcnt 
quadratic form H : T*M — > R - the system's kinetic energy P(p) = T(p) = i < p,p >, in local 
canonical coordinates x t ,pi € U p , U p open in T*M, given by 

H{P) = \g i0 {x,m)pip , (7) 

where g l i(x,rn) — g^ (x,m) denotes the inverse (contravariant) material metric tensor 

v— -\ Ox* &x^ 

g v {x ,m)=^m x 6 rs — — . 

x=i 

T*M is an orientable manifold, admitting the standard volume form 

, W(W + 1) 

Wff ~ AH 

For Hamiltonian vector-field, Xh on M, there is a base integral curve 7o(t) = Pi(t)) iff 

7o(t) is a geodesic, given by the one-form /orce equation 

Pi=Pi+T i jk gi l g km p lPm = 0, with x k = g ki Pi . (8) 

The l.h.s pi of the covariant momentum equation ([5]) represents the intrinsic or Bianchi covariant 
derivative of the momentum with respect to time t. Basic relation pi — defines the parallel 
transport on T , the simplest form of human-motion dynamics. In that case Hamiltonian vector- 
field Xh is called the geodesic spray and its phase-flow is called the geodesic flow. 



for any two forms a, (3 £ f2 p (Jf ). 

7 The Lie derivative ZuOi of p— form a along a vector-field u is defined by Cartan's 'magic' formula (see [4][5]): 

SmCt = u\dct + d(u\a). 

It satisfies the Leibnitz relation 

£ u {a A 0) = to A/) + ttA £u/3. 

Here, the contraction J of a vector-field u = u^d^ and a p— form a = axi...\ p da: 1 A • • ■ A dx x p on a biomechanical 
manifold X is given in local coordinates on X by 

u\a = u'*a ll x 1 ...x p _ 1 dx Xl A ■ ■ ■ A dx x P- 1 . 

It satisfies the following relation 

u\{aA(3) = «JaA/3 + (-l) |a| QAuJ/3. 
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For Earthly dynamics in the gravitational potential field V : M — > R, the Hamiltonian H : 
T*M ->K0 extends into potential form 

with Hamiltonian vector-field -X"h = JVH still defined by canonical equations ([6]). 
A general form of a driven, non-conservative Hamiltonian equations reads: 

x^d^H, Pi =Fi-d x iH, (9) 

where Fi = -Fi(i, x,p) represent any kind of joint-driving covariant torques, including active neuro- 
muscular-like controls, as functions of time, angles and momenta, as well as passive dissipative and 
elastic joint torques. In the covariant momentum formulation l[8"]). the non-conservative Hamilto- 
nian equations ([9]) become 

Pi = Pi + T) k g 3l g km pip m = Fi, with x k = g kl Pi . 

The general form of autonomous Hamiltonian biomechanics is given by dissipative, driven 
Hamiltonian equations on T*M: 

x = -^T + ^Ti ( 10 ) 
o P i o Pi 

^ dH dR 

x\0) = 4, Pi(0)=p°i, (12) 

including contravariant equation ([T0|) - the velocity vector-field, and covariant equation (p~T|) - 
the /orce 1-form (field), together with initial joint angles and momenta (I12j) . Here i? = R{x,p) 
denotes the Raileigh nonlinear (biquadratic) dissipation function, and Fi = Fi(t, x,p) are covariant 
driving torques of equivalent muscular actuators, resembling muscular excitation and contraction 
dynamics in rotational form. The velocity vector-field (TIT))) and the force 1— form Ijlip together 
define the generalized Hamiltonian vector-field Xh\ the Hamiltonian energy function H = H(x,p) 
is its generating function. 

As a Lie group, the biomechanical configuration manifold M = FT • SE(3y is HausdorffH There- 
fore, for x = (x l , p{) G Up, where U p is an open coordinate chart in T*M, there exists a unique 
one-parameter group of diffeomorphisms cf> t : T*M — > T* M , that is the autonomous Hamiltonian 
phase-flow: 

^ : T*M -,T*M :(p(0),x(0))^(p(t),x(t)), (13) 
{4>t ° <f>s = 4>t+si 4>q = identity), 

given by (fTTM2|) such that 

j t | t=0 <j> t x = JVH(x). 



8 That is, for every pair of points Xx,%2 G M, there are disjoint open subsets (charts) XJ\,XJi C M such that 
x\ £ Ui and X2 £ C/2- 
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4 Time— Dependent Hamiltonian Biomechanics 



In this section we develop time-dependent Hamiltonian biomechanics. For this, we first need to 
extend our autonomous Hamiltonian machinery, using the general concepts of bundles, jets and 
connections. 

4.1 Biomechanical Bundles and Jets 

While standard autonomous Lagrangian biomechanics is developed on the configuration manifold 
X, the time-dependent biomechanics necessarily includes also the real time axis R, so we have 
an extended configuration manifold Ixl. Slightly more generally, the fundamental geometrical 
structure is the so-called configuration bundle n : X — > R. Time-dependent biomechanics is thus 
formally developed either on the extended configuration manifold K x X, or on the configuration 
bundle 7r : X — > R, using the concept of jets, which are based on the idea of higher-order tangency, 
or higher-order contact, at some designated point (i.e., certain anatomical joint) on a biomechanical 
configuration manifold X. 

In general, tangent and cotangent bundles, TM and T*M, of a smooth manifold M, are special 
cases of a more general geometrical object called fibre bundle, denoted ir : Y — > X, where the word 
fiber V of a map tt : Y — > X is the preimage tt~ 1 (x) of an element x G X. It is a space which 
locally looks like a product of two spaces (similarly as a manifold locally looks like Euclidean space), 
but may possess a different global structure. To get a visual intuition behind this fundamental 
geometrical concept, we can say that a fibre bundle Y is a homeomorphic generalization of a 
product space X x V (see Figure [2]), where X and V are called the base and the fibre, respectively. 
tt : Y — > X is called the projection, Y x = n^ 1 (x) denotes a fibre over a point x of the base X, while 
the map / = ir^ 1 : X — > Y defines the cross-section, producing the graph (x, f(x)) in the bundle 
Y (e.g., in case of a tangent bundle, / = x represents a velocity vector-field). 



bundle Y^XxV 




Figure 2: A sketch of a locally trivial fibre bundle Y « X x V as a generalization of a product 
space I x V; left - main components; right - a few details (see text for explanation). 

More generally, a biomechanical configuration bundle, ir : Y — ► X, is a locally trivial fibred 
(or, projection) manifold over the base X. It is endowed with an atlas of fibred bundle coordinates 
(x x ,y l ), where (x A ) are coordinates of X. 

Now, a pair of smooth manifold maps, f\, f^ ■ M — > N (see Figure[3]), are said to be k— tangent 
(or tangent of order k, or have a fcth order contact) at a point a; on a domain manifold M , denoted 
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by h ~ /a, iff 



/i(x) = /2(^) called — tangent, 
d x fi(x) = d x f 2 (x), called 1 - tangent, 
d X xfi(x) = d xx f 2 (x), called 2 - tangent, 
etc. to the order fc 

In this way defined k—tangency is an equivalence relation. 




i n 

Figure 3: An intuitive geometrical picture behind the A; —jet concept, based on the idea of a higher- 
order tangency (or, higher-order contact). 

A k—jet (or, a jet of order k), denoted by j*f, of a smooth map / : M — > 2V at a point a; G M 
(see Figure[3]), is defined as an equivalence class of k— tangent maps at x, 

j*f : M -> N = {/' : /' is k - tangent to / at a:}. 

For example, consider a simple function f : X — > Y, x ^ y = f(x), mapping the X— axis into 
the Y — axis in R 2 . At a chosen point x S X we have: 
a 0— jet is a graph: (x, f(x)); 
a 1— jet is a triple: (x, f(x), f'(x)); 
a 2— jet is a quadruple: (a;, /(x), /'(#), f"(x)), 
and so on, up to the order fc (where /'(x) = etc). 

The set of all k— jets from j*/ : X — » F is called the A; —jet manifold J k (X, Y). 

Formally, given a biomechanical bundle Y — > X, its first-order je£ manifold J X Y comprises 
the set of equivalence classes j^s, x 6 X, of sections s : X — > F so that sections s and s' belong 
to the same class iff 

Ts \t so x= Ts' | Tx x ■ 

Intuitively, sections s,s' € j^s are identified by their values s l (x) — s' l (x) and the values of their 
partial derivatives d fl s t (x) — d^s' l (x) at the point x of X. There are the natural fibrations [15j 

7Ti : J X Y 3 jls i — >i£ X 1 7T i : J X Y 3 j x s h-> s(x) G F. 

Given bundle coordinates (x x ,y l ) of F, the associated jet manifold J F is endowed with the 
adapted coordinates 

(x\y\yl), (y\y\)(jls) = (s*(x),d x s*(x)), = ( dfl + y^y'*. 
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In particular, given the biomechanical configuration bundle M — > R over the time axis R, 
the 1—jet space J 1 (R, M) is the set of equivalence classes j^s of sections s 1 : K — > M of the 
configuration bundle Af — > R, which are identified by their values s l (i), as well as by the values 
of their partial derivatives dts z = dts 1 ^) at time points t g K. The 1—jet manifold J 1 (R, M) is 
coordinated by (t, x l ,x % ^, that is by (time, coordinates and velocities) at every active human joint, 
so the 1-jets are local joint coordinate maps 

jjs-.R^M, t^(t,x l ,x l ). 

The second-order jet manifold J 2 Y of a bundle Y — > X is the subbundle of J 2 Y — ► J X Y 
defined by the coordinate conditions y l Xfl — y^ A . It has the local coordinates (x x ,y l ,y\,y x<fJ ^) 
together with the transition functions [To] 

dx a 

y% = g^i 9 * + y* d J + yia 9 j)y'x- 

The second-order jet manifold J 2 Y of Y comprises the equivalence classes j 2 s of sections s of 
Y — > X such that 

1/1(Jx«) - dxs l (x), yl^s) = d^(x). 

In other words, two sections s, s' € j 2 s are identified by their values and the values of their first 
and second-order derivatives at the point x G X. 

In particular, given the biomechanical configuration bundle M — > R over the time axis R, 
the 2— jet space J 2 (R, M) is the set of equivalence classes jfs of sections s 1 : R — > M of the 
configuration bundle 7r : M — > R, which are identified by their values s l (i), as well as the values 
of their first and second partial derivatives, dts 1 = dts l (t) and dus 1 = dtts l (t), respectively, at 
time points t 6 R. The 2— jet manifold J 2 (R, M) is coordinated by (t, x l , x 1 , £*), that is by (time, 
coordinates, velocities and accelerations) at every active human joint, so the 2-jets are local joint 
coordinate map^l 

j 2 s:R-^M, 1 1-» (t,x\x\x l ). 

4.2 Nonautonomous Dissipative Hamiltonian Dynamics 

We can now formulate the time-dependent biomechanics in which the biomechanical phase space 
is the Legendre manifold II, endowed with the holonomic coordinates (t,y l ,pi) with the transition 
functions 

II admits the canonical form A given by 

A = dpi A dy l A dt <g> d t . 

We say that a connection 

T = dt ® (d t + fd t + lt d l ) 

9 For more technical details on jet spaces with their physical applications, see |15l I16| ). 
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on the bundle II — ► X is locally Hamiltonian if the exterior form 7J A is closed and Hamiltonian if 
the form 7JA is exact [15] , A connection 7 is locally Hamiltonian iff it obeys the conditions: 

SV - «9V = 0, da - = 0, drf + 0*7j = 0. 

Note that every connection Y = dt ® (dt + r z di) on the bundle Y — > X gives rise to the 
Hamiltonian connection T on n — >1, given by 

f = dt ® (dt + T% - OA ',,,<)■■ ). 

The corresponding Hamiltonian form Hr is given by 

H T = pidy 1 - piVdt. 

Let H be a dissipative Hamiltonian form on n, which reads: 

H = p t dy l - Hdt = p t dy l - p % Y l dt - H T dt. (14) 

We call TL and TL in the decomposition (|14[) the Hamiltonian and the Hamiltonian function re- 
spectively. Let 7 be a Hamiltonian connection on II — > X associated with the Hamiltonian form 
HID). It satisfies the relations QHQI] 

7J A = dpi A dy l + 7^7/ A dt - -y l dpi A dt = dH, 

7* = d l U, lt = -diH. (15) 

From equations (fTS)) we see that, in the case of biomechanics, one and only one Hamiltonian 
connection is associated with a given Hamiltonian form. 

Every connection 7 on n — > X yields the system of first-order differential equations: 

f = l\ Pi = 7i- (16) 

They are called the evolution equations. If 7 is a Hamiltonian connection associated with the 
Hamiltonian form H (fill , the evolution equations (|TB)) become the dissipative time- dependent 
Hamiltonian equations: 

V* = d l H, Pi = -diH. (17) 

In addition, given any scalar function / on n, we have the dissipative Hamiltonian evolution 
equation 

d H f = (d t + d'Hdi - dm 1 ) f, (18) 

relative to the Hamiltonian 7i. On solutions s of the Hamiltonian equations (|17j) . the evolution 
equation (fl8|) is equal to the total time derivative of the function /: 

s*d H f=^-(fo S ). 
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4.3 Time— Dependent Biomechanics 

The dissipative Hamiltonian system (fT"7l) - (fT8f is the basis for our time & fitness-dependent biome- 
chanics. The scalar function / in (fT5)) on the biomechanical Legendre phase-space manifold II is now 
interpreted as an individual neuro-muscular fitness function. This fitness function is a 'determi- 
nant' for the performance of muscular drives for the driven, dissipative Hamiltonian biomechanics. 
These muscular drives, for all active DOF, are given by time & fitness-dependent Pfaffian form: 
Fi = Fi(t,y,p, /). In this way, we obtain our final model for time & fitness-dependent Hamiltonian 
biomechanics: 

if = &H, 
pi = Fi - diH, 
d H f = (8 t + d l Hd i -dSHd i )f. 

Physiologically, the active muscular drives Fi = Fi(t,y,p, /) consist of [H[3]) : 

1. Synovial joint mechanics, giving the first stabilizing effect to the conservative skeleton 
dynamics, is described by the (y,y)-form of the Rayleigh-Van der Pol's dissipation function 

n 

i=l 

where on and (i i denote dissipation parameters. Its partial derivatives give rise to the viscous- 
damping torques and forces in the joints 

P° int = dR/dy\ 

which arc linear in y l and quadratic in y % . 

2. Muscular mechanics, giving the driving torques and forces F™ usc = Tl ausc {t,y,y) with 
(i = 1, . . . , n) for human biomechanics, describes the internal excitation and contraction dynamics 
of equivalent muscular actuators |12j . 

(a) The excitation dynamics can be described by an impulse force-time relation 
pvra-p = pQ(i_ e -t/r^ if stimulation > 
F un P = p o e -t/r, if stimulation =0, 

where F® denote the maximal isometric muscular torques and forces, while Tj denote the associated 
time characteristics of particular muscular actuators. This relation represents a solution of the 
Wilkie's muscular active-state element equation [13] 

fi + Tfi = TSA, (i(0) = 0, 0<5<1, 

where /i = /j,(t) represents the active state of the muscle, T denotes the element gain, A corresponds 
to the maximum tension the element can develop, and S = S(r) is the 'desired' active state as a 
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function of the motor unit stimulus rate r. This is the basis for biomechanical force controller. 

(b) The contraction dynamics has classically been described by the Hill's hyperbolic force- 
velocity relation [14] 

pHUi = ( F ? b * ~ 'V' -•'>') 
(Siji* + h) ' 

where a, and hi denote the Hill's parameters, corresponding to the energy dissipated during the 
contraction and the phosphagenic energy conversion rate, respectively, while 5ij is the Kronecker's 
S— tensor. 

In this way, human biomechanics describes the excitation/ contract ion dynamics for the ith 
equivalent muscle-joint actuator, using the simple impulse-hyperbolic product relation 

T^ usc {t,y,y) = Fi mp x F t mu . 

Now, for the purpose of biomedical engineering and rehabilitation, human biomechanics has 
developed the so-called hybrid rotational actuator. It includes, along with muscular and viscous 
forces, the D.C. motor drives, as used in robotics 

H 0b0 = *fc(0 - JkVk(t) - B k y k (t), with 
Wfc(f) + Rkik(t) + ChVk{t) = Uk{t), 

where k — 1, . . . , n, ifc(t) and Uk(t) denote currents and voltages in the rotors of the drives, Rk, Ik 
and Cfe are resistances, inductances and capacitances in the rotors, respectively, while Jk and Bk 
correspond to inertia moments and viscous dampings of the drives, respectively. 

Finally, to make the model more realistic, we need to add some stochastic torques and forces: 

Ff och = B. l] {y*(t),t]dW>(t), 

where Bij [y(t) , t] represents continuous stochastic diffusion fluctuations, and W 3 (t) is an N— variable 
Wiener process (i.e., generalized Brownian motion) [6], with 

dW j (t) = W 3 {t + dt) - W 3 (t), (for j = 1, . . . ,n = no. of active DOF). 

5 Conclusion 

In this paper we have proposed the time-dependent Hamiltonian form of human biomechanics. 
Starting with the Covariant Force Law: Fi — niija J on the biomechanical configuration manifold 
M, we have first developed the autonomous Hamiltonian biomechanics: 

±i _dH OR . _ p dH + dR 

dpi dpi ' ^ % 1 dx l dx 1 ' 

on the symplectic phase space that is the cotangent bundle TM of M. Then we have intro- 
duced powerful geometrical machinery consisting of fibre bundles and jet manifolds associated to 
the biomechanical manifold M. Using the jet formalism, we derived time-dependent, dissipative, 
Hamiltonian equations: 

f = shi, p.- = fi - cw, 
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together with the fitness evolution equation: 

d H f = (d t + d'Hdi - d.HO 1 ) /. 
for the general time-dependent human biomechanical system. 
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